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Abstract 

Renormalization is cast in the form of a Lie algebra of infinite triangular matrices. By expo- 
nentiation, these matrices generate counterterms for Feynman diagrams with sub divergences. 
As representations of an insertion operator, the matrices are related to the Connes-Kreimer 
Lie algebra. In fact, the right-symmetric nonassociative algebra of the Connes-Kreimer in- 
sertion product is equivalent to an "Ihara bracket" in the matrix Lie algebra. We check our 
results in a three-loop example in scalar field theory. Apart from possible applications in 
high-precision phenomenology, we give a few ideas about possible applications in noncom- 
mutative geometry and functional integration. 
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1 Introduction 



It has been known for several years now that renormahzation is governed by certain algebras 
known as Hopf algebras (for a review, see Kreimer IHD])- Unfortunately, since Hopf algebra is 
largely unfamiliar to physicists^, this intriguing fact has not been widely appreciated. Here, 
we make an attempt to help remedy this situation, by focusing instead on an equivalent 
Lie algebra. The existence of such a Lie algebra is guaranteed by the famous Milnor-Moore 
theorem (see e.g. [Ej), a fact that was exploited by Connes and Kreimer already in 
This theorem says, loosely speaking, that any graded Hopf algebra can be written in terms of 
a Lie algebra. The mere fact that Lie algebra is much more familiar to physicists than Hopf 
algebra makes this direction useful to explore. In addition, there are certain computational 
advantages that may prove important; the Lie algebra can be expressed in terms of (in general 
infinite) matrices, so renormahzation becomes a matter of matrix algebra. The matrix Lie 
algebra is the subject of this work. We explore the application to renormahzation and give 
a number of examples. 

Via the homomorphism between the Hopf algebra of renormahzation and the Hopf 
algebra of coordinates on the group of formal diffeomorphisms in noncommutative geometry 
Pm 122] , we also expect to be able to apply the matrix framework in noncommutative geom- 
etry terms, but this direction is postponed to future work. (Reasons why noncommutative 
geometry could be relevant to physics abound jHl El IHH EE]-) There is also another, more 
speculative way in which this matrix Lie algebra may be applied, explained in the conclu- 
sion, section El First, let us step back and review some existing directions in the literature. 
Partly through indications of recommended references, we will attempt to maintain some 
accessibility for mathematicians who are not completely familiar with physics jargon. 

2 Background 

Perturbative renormahzation was regarded by most physicists as finished with a theorem by 
Bogoliubov, Parasiuk, Hepp, and Zimmermann (the BPHZ theorem), refined by Zimmer- 
mann in 1970. Original references for this theorem are |2], and a clear textbook version is 
given in [16^. Already hidden within the BPHZ theorem, however, was algebraic structure 
belonging to branches of mathematics that were probably largely unknown to physicists in 

■^except in completely different contexts, as symmetries of some nonlinear sigma models in low-dimensional 
supergravity. A basic Hopf algebra reference is • 
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the 1970s. The Hopf algebra of renormahzation was first displayed by Kreimer in 1997, 
partially as a result of his excursions into knot theory PHJ. In 1999, Broadhurst and Kreimer 
performed a Hopf-algebraic show of strength by computing contributions to anomalous di- 
mensions in Yukawa theory to 30 loops |4|'^. 

On a practical level, phenomenologists have been opening their eyes to Hopf algebra as 
a potential time-saver and organizing principle in massive computations, such as the 5-loop 
computation of the beta function in quantum chromodynamics (QCD). Such high-precision 
analytical calculations are needed for comparisons to numerical calulations in lattice QCD, 
since the energy scales accessible in these lattice computations (and in many experiments) 
are often so low that perturbative QCD is only starting to become a valid approximation 
(see e.g. ^1 E]). Even when one is not seeking to push calculations to such heights of 
precision, savings in computer time through better algorithms could be welcome for routine 
calculations as well. 

On the more mathematical side, a series of articles by Connes and Kreimer 1201 
took the Hopf algebra of renormahzation through some formal developments, amongst other 
the identification of the Lie algebra of the dual of the aforementioned Hopf algebra, through 
the Milnor-Moore theorem (see also fO] for related mathematical developments). This has 
later been phrased as a study of operads^, which we will not go into here, but which can 
be seen as complementary to our discussion. Very recently, some geometric aspects were 
explored in [T^ . 

We believe that representations of the dual Lie algebra can be useful in their own 
right, so this work is about one such representation. The Lie bracket is simply given by 
inserting one graph into another, and subtracting the opposite insertion. The process of 
insertion may be represented by the multiplication of (in general infinite) matrices, hence 
matrix Lie algebra. 

Although we give examples in scalar field theory, the algebra representation on graphs 
is independent of the specific nature of the action, except that we assume the interactions 
are quartic. We could equally well consider fermions, vector fields, or higher-spin fields 
with quartic interactions. But surely Ward identities, ghosts and so on make gauge theory 
computations very different from those of scalar field theory? This question can be answered 
on two levels. First, the distinction between Feynman graphs and their values. In performing 
^This was later improved further 5 . 

^For instance, a simple example of an operad is a collection of maps from all tensor products of an algebra 
A into A itself, such as A ® A ® A ^ A, with some compatibility requirements. See e.g. for a precise 
definition. 
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field theory computations, one routinely identifies graphs and their values. In the Connes- 
Kreimer approach one explicitly separates symmetries of the algebra of Feynman graphs 
from symmetries of Feynman integrals; the integrals are thought of as elements of the dual 
of the space of graphs. The dual has symmetries of its own, "quasi-shufiles" [22] , obviously 
related but not identical to symmetries on graphs. The algebraic approach can help with 
computing Feynman integrals, but its greatest strengths so far has been in the organization 
of counterterms, providing algebraic relations between seemingly unrelated quantities. The 
organization of counterterms is an essentially graph-theoretic problem, since the hierarchy 
of subdivergences can be read off from the topology of a graph without knowing precisely 
what fields are involved. In the extreme case of "rainbow" or "ladder" diagrams only, the 
Feynman integrals themselves are easily iterated to any loop order (as in |4), but the BPHZ 
subtractions quickly become a combinatorial mess if the algebraic structure is ignored. 

The attitude of separating the combinatorial and analytical problems (space of graphs 
and its dual) is also useful for computer implementation. Let us take one example: the 
Mathematica package FeynArts j21I, which can automatically compute amplitudes up to 
one loop in any renormalizable theory, and using additional software, up to two loops [25 . 
This particular software package calculates amplitudes by first writing down scalar graphs 
with appropriate vertices, then generates counterterm graphs, then at a later stage lets the 
user specify what fields are involved, either gauge fields, spinors or scalars. Since the number 
of graphs grows very quickly with loop order, it will become important to take advantage 
of algebraic structure if calculations at very high precision are to be feasible in a reasonable 
amount of computer time. 

Apart from the direct application of the matrix Lie algebra to precision computations, 
there are other possible indirect applications. Such applications may be found in noncom- 
mutative geometry through the result of [201. A currently more speculative application is 
to functional integration. These directions are not explored in detail in this work, but see 
section [TUl 

3 Graph Summary 

We consider 0*^ theory in four spacetime dimensions for ease of exposition (see the previous 
section for comments on adaption to other theories). In 0^ theory in four dimensions, we 
consider graphs with a quartic interaction vertex^ and two and four external legs (see section 
^i.e. the interaction vertex is four-valent. 
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Table 1: Graphs with L loops and E external legs. 



Elfor a reminder of why we do not need e.g. six external legs in this theory in four dimensions). 
All such graphs with up to three loops are summarized in table In this table, "crossing" 
refers to graphs related to previously given graphs by crossing symmetry, such as the two 
graphs in parenthesis for L = 1, E = 4; they are simply found by turning the first graph on 
its side (second graph), and then crossing the two lower external legs (third graph). ^ On a 
computer, graphs can be stored as lists. One example of such lists, adapted to the context 
of Feynman graphs, is given at the end of section ^ 

4 The Grafting Operator 

The Lie algebra of graphs comes from inserting (grafting) one Feynman graph into another. 
First a remark about terminology. In standard Feynman graphs, each internal line represents 
a propagator. We will only work with amputated graphs, where the external legs do not 
represent propagators. Now, in a theory with quartic interactions, one can consider inserting 
graphs with either two external legs (into an internal line) or four external legs (into an 
interaction vertex). Thus, the middle of an internal line and the interaction vertex itself are 
eligible as insertion points. If desired, we could represent each insertion point by a dot, as 
^ These graphs were drawn using FeynMF I5gj . 

®The reader unfamiHar with Feynman diagrams may want to note that the third graph in L = 1, = 4 
still has only two vertices: the two lines crossing below the loop are going one above the other. 
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in the graph on the left in fig. Q As in standard usage — the graph on the right in fig. ^ — we 
leave those dots implicit. We alert the reader that the authors of JHl call the dots on lines 




Figure 1: Insertion points are implicit. 

"two-point vertices" . To avoid confusion with some of the physics terminology, we prefer to 
talk about insertion points instead of two-point vertices. 

Here is a precise definition of the insertion, or grafting, operator. The practically- 
minded reader can refer to the insertion tables in appendixElto get a quick idea of the action 
of this operator. Enumerate the insertion points of a graph t by Pi{t), the total number of 
insertion points of t by N{t) and their valences (number of "prongs", i.e. lines sticking out 
from the vertex) by v{pi{t)). Here, the valence is always either 2 or 4. If ti and t2 are graphs, 
and the number f (^2) of external legs of ^2 is equal to the valence v{pi{ti)) of the insertion 
point Pi{ti) of ti, then we define s^^''*^^ to be insertion of the graph t2 at Pi{ti) by summing 
over all permutations of the v{pi{ti)) propagators to be joined. Then the grafting operator 
■5*2 is 

. Niti) 

The total number of graphs A/" (including multiplicities of the same graph) created by inser- 
tion of t2 into ti is 

N{ti) 

= ^HMii))MMt2)MPi{ti)y- 
1=1 

We call ti the "object graph" and t2 the "insertion graph" . Often, we will use parenthesis 
to delineate the insertion graph: s(t2)^i- Finally, we define s{t) to be linear: 

s{ati + bt2) = as{ti) + bsit2) a,b E C . 

A remark about normalization: the normalization in is a mean over insertions. It gives 
Sjl = t, where 1 represents tree- level graphs, with either 2 or 4 external legs. (When one 
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needs to make a distinction, one can use I2 and I4.) On the other hand, we have sit 7^ t, in 
fact 

Sit = N{t)t = {3L{t) ±l)t (tj^l) 

where L{t) is the number of loops of t and the upper (lower) sign refers to graphs with 4 (2) 
external legs, respectively. 

We can assign a grading (degree) to the algebra by loop number L, since then obvi- 
ously 

deg st^ti = degta + degti , 

so within the set of graphs with a fixed number of external legs, Si acts as our grading 
operator. Being the only st which is diagonal, we will exclude si from the algebra, but it 
can be adjoined whenever needed as in pUj . 

As an alternative to this construction of the grafting operator, one could consider a 
restricted grafting operator, where permutations of the legs of the insertion which lead to 
the same topologies are represented by separate operators. In this work, we concentrate on 
exploring the full operator. 

4.1 Nested divergences 

It is of particular interest to know when a graph has nested divergences, i.e. when two 
divergent loops share the same propagator. For any graph, two lines are said to be indis- 
tinguishable if they are connected to the same vertex, and distinguishable if they are not. 
Whenever we insert a (loop) graph with four external legs at a certain vertex and two dis- 
tinguishable external prongs of the insertion graph are connected to indistinguishable legs 
of another vertex, we create nested divergences. Of course, any nested divergences already 
present in the insertion graph are preserved. 

Here is an example. Consider the insertion in fig. [21 It is an insertion of ^ into the 
upper vertex of that same graph. With the above use of language, 5 and 6 are equivalent, 
and so are 9 and 10, but not e.g. 5 and 7. When we join 3-5, 4-6, 7-9, 8-10, we create 
which only has two disjoint divergent loops. But when we join 3-5, 4-8, 7-9, 6-10, we create 
which has nested divergences. 

As a side remark, figure El also lets us recall how X is easily represented in a format 
amenable to computer processing. The bottom vertex can be stored as [9,10,11,12], and 
the propagator going between 5 and 7 carrying loop momentum p can be stored as [5, 7,p]. 
This way, the list of vertices contains three sublists, and the list of propagators contains four 



7 



1 2 




11 12 



Figure 2: The creation of overlapping divergences 

sublists, which together completely specify the graph. This is also how graphs are stored 
in FeynArts |21]. Of course, there are many other equivalent representations, for instance 
using the relation between rooted trees and parenthesized words j^]. 



5 Lie Algebra 

We want to show that the grafting operators St form a Lie algebra C To this end, we first 
show a certain operator identity. Indeed for graphs ti, t2, define 

[tut2] = s{h)t2-s{t2)h (2) 

then we shall show that the commutator of the operators s(ti), 5(^2) is equal to s{[ti,t2]): 

[S{h),s{t2)]=s{[ti,t2]) , (3) 

as an operator identity acting on graphs. This identity is analogous to the Ihara bracket of 
the Magnus group, familiar from number theory (see e.g. j35j). 

Here is a sketch of the proof of 0- Consider two insertions into a graph ^3, s^f^^^ and 
^Pj{*3)^ before summing over insertion points Piit^) and Pj{t^). Let us separate insertions into 
"mutually local" {lipiit^) = Pjit^)) and "mutually nonlocal" (otherwise). This is illustrated 
in figure El The "mutually nonlocal" insertions clearly commute, since they modify different 
points, so these contributions to St-^ and cancel in the bracket. We are left with the 
"local" insertions, i.e. 

Stit2 - Stji =: [ti,t2] 
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Figure 3: Difference between mutually nonlocal (a) and mutually local insertions (b). 



This concludes the outline of the proof. Using the correspondence in section 15.11 a detailed 
proof will appear in p. 

We remark that the identity trivially holds acting on 1: 

[S{h),s{t2)]l = S{ti)s{t2)l - S{t2)s{ti)l = S{h)t2 - ^(ta)^! = [^1,^2] • 

It is easy to check that the bracket [^1,^2] = s{ti)t2 — s{t2)ti satisfies the Jacobi 
identity. By writing down three copies of the Ihara bracket identity and adding them we 
find by linearity 

s{[tut2])t3 + s{h)[t2,ti] + cyclic = 

or 

[tl,[t2,h]]+CJclic = . 

The Lie algebra is graded by loop order: C = C^^^ © C^'^^ © C'^^^ © . . .. and 

J^(m)^^(„)j C . 

Let us see what the Lie bracket looks like explicitly in our field theory with four- 
vertices. To make the notation more economic, we suppress any graph which is related by 
crossing symmetry to one that is already included. That is, since ^, )o( and (see table H]) 
all represent the same function of the different Mandelstam variables^ s, t, and u, we will 
only display ^ out of those three diagrams. We can always use a symmetric renormalization 

^Here s — (pi +^2)^, t= {p'l ~ pi)^ and u — — pi)^, where pi, p2 are incoming momenta, andp^, p'2 
are outgoing momenta. 
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point, i.e. we can define the theory at some mass scale M, meaning s = t = u = -M^ at this 
scale, in which case counterterms are the same for all three diagrams mentioned previously. 

At the one-loop level, we have the following diagrams: — , X) ^i^d ^. Consulting 
the insertion tables in the appendix, or equivalently the matrix representation which we 
introduce below, we find the commutators of tree-level diagrams with one-loop diagrams: 

[£(0)^£{l)] . 

[-,-2-] = si-)lL - .s{SL)- 

= 2SL-SL 

= IL 

[x,x] = s(xm-s{^)x 

The only nonvanishing commutator of two one-loop graphs is 

c £(2) : 

= ^ 0{h). (4) 

where by — ^ we mean evaluation of the graphs. Here we have restored h to connect with 
the semiclassical approximation; an L-loop graph is of order h^~^. At this point we note that 
the combinatorial factors 1/3, 2/3 and so on are not the usual symmetry factors of graphs, 
i.e. the multiplicity of operator contractions generating the same graph. 

Including two-loop graphs in the algebra, the following graphs come into play: X) S!) 
^ and We can now consider a commutator of a one-loop and a two-loop graph, which 
creates a sequence of three-loop graphs: 



[£l",£W] c£« : 



5X + |S + |S 



3 
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which may be approximated by zero if we only keep graphs up to order h. Thus, to this 
order, all the "new" graphs with two loops are central elements of the algebra. 

In general, if we consider L loops (order h^^^), the commutator at order is negligi- 
ble, so only commutators of graphs of loop orders Li, L2 where L1 + L2 < L are nonvanishing. 
The effect of going from L to L + 1 is then to add central graphs as new generators, which 
appear in lower-order commutators with one extra factor of h. In other words, the new 
graphs deform the L-loop algebra by h. 

Now that we have displayed a few commutators, here are two examples of how the 
Ihara bracket identity works. Let us first consider a one- loop/one- loop commutator, acting 
on the bare propagator. 

s(s(n).SL)--s{sWS.)- 
1 2 

--a_-F --e- 

3 3 

[s(X),s(iL)]- = sC^)s{SL)--s(SL)sC^)- 
1 2 

= -_a_-F--e-. 

3 3 



whereas 



At this level, the identity is quite trivial. Let us therefore also check the relation at the 
three-loop level: 

ls(5),s(-a-)lS = s(5)s(j!.) 5 - 5 

= s(S)(2$) - (5 X + 5 g + 1 
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but 

s([5,-a-])S = s('i-a-+|^-2g') 5 

Here we see that cancellation of graphs ^and ^was necessary in the first bracket, which 
is to be expected since these two graphs are generated by "mutually nonlocal" insertions, as 
defined above. 

To summarize, if we know the commutator of two insertions, the identity (jS)) gives us 
the insertion of the commutator. A more concrete expression of this will be seen in section 

Going back to the question of normalization, not any combination of normalizations 
will preserve the Ihara bracket, of course, as can easily be seen from the example (jH)). One 
way to preserve it is to simultaneously drop overall normalization and multiplicities, but 
without multiplicities the grafting operator loses some of its combinatorial information (it 
ceases to count the number of ways a graph can be generated). Although from a physics point 
of view there is no freedom in the choice of normalization in the total sum of renormalized 
graphs — for given Feynman rules — it can be useful to consider different normalizations in 
these intermediate expressions. 

5.1 Connection with the Star Operation 

The Ihara bracket identity can be compared to the star product discussed by Kreimer 
Here is the connection between our operator St and Kreimer's star insertion -k: 

ti-kt2 ^ S{t2) ti 

where the ~ alerts the reader that we use a different normalization, (Q). 
Furthermore, 

{tlkt2)kh ~ S{h)s{t2)ti (6) 
tl^{t2^h) ~ S{s{h)t2)ti. (7) 
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The coproduct A of the Connes-Kreimer Hopf algebra is coassociative, so the dual 
Hopf algebra should have an associative multiplication 0. One would then naively expect 
the dual Lie algebra to have an associative multiplication as well. However, comparing (jHI) 
and dZj), the star operation is not associative. This is because the original Hopf algebra con- 
tains disconnected graphs, but only connected graphs are generated by the star product (or 
equivalently, the grafting operator St). When the multiplication is truncated to connected 
graphs, it is no longer associative. 

The deviation from associativity can be measured by the associativity defect: 

a{ti I t2, ts) ^h-k {t2 -k ts) - {ti -k t2) -k tg . (8) 

When this is nonzero, as in any nonassociative algebra, it is interesting to consider the 
"right-symmetric" algebra for which 

a{h I ^2,^3) = I ^3,^2) (9) 

(a "left-symmetric" algebra can equivalently be considered). This leads naturally to an 
algebra of the star product [201! which is also common usage in differential geometry. Writing 
the right-symmetric condition (jH)) explicitly in terms of the associativity defect (jH)), the 
algebra is specified by 

ti -k {t2 -k ts) - {ti -k ta) 'kts = ti'k (tg -k ta) - (tl -k tg) -k . (10) 

This is known as a pre-Lie (or Vinberg) algebra. A pre-Lie algebra yields a bracket [^1,^2] = 
ti-kt2—t2kti that automatically satisfies the Jacobi identity, hence "pre-Lie". In our notation, 
using eqs. ©, ©j the pre-Lie algebra identity (fTUj) can be rewritten as the operator identity 
(acting on arbitrary tg): 

[S(tl),s(t2)] =s([ti,t2]) , (11) 

that we have already shown directly for the grafting operator st- Thus the two descriptions 
are equivalent, as expected. 

5.2 Matrix Representations 

Graphs are graded by their number of loops, but there is no canonical ordering within each 
graded subspace. Given some ordering of graphs, the grafting operator Sj can be represented 
as a matrix. This matrix will be indexed by graphs. Now, graphs are inconvenient to use as 
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subscripts, so we use a bra-ket notation instead: represent s(ti)t2 by 

Sih)\t2) = Y,\t3){t3\s{tl)\t2) 

where all graphs are to be thought of as orthogonal: (ti|t2) = for ti ^ t2- 
An example: the insertion tables (consulting appendix El gives 

which is to be compared to the expansion in a complete set 

s{yi)\^) = ... + {x\sm^)\x) + msm^)\^) + Msm^)\^) 

so we read off 

( X \sm X ) = 2/3 , \sm X ) = 2/3 , \sm x ) = 2/3 . 

giving the second column of the matrix: 



/ 0\ 


/o\ 






1/3 




1 






2/30000 









2/3 


2/30000 









2/3 


2/30000 









2/3 


V 0/ 


\o/ 




V J 



,2 2 2 

= 3X+3X+3X (12) 



where we have ordered the graphs as X)5!? X) S^jSk, and = denotes "represents". In 
the appendix we give first insertion tables up to three- loop order (appendix^, then a few 
examples of matrices that represent the grafting operators (appendix EI). These matrices 
are easily extracted from the insertion tables, and are indexed by graphs. 

A word about practical implementation: because of our space-saving convention of 
suppressing graphs related by crossing symmetry, each number in the above 6x6 matrix is 
a 3 X 3 unit matrix. The exception is the first column, which has as second index the single 
tree-level diagram and so is a 3 x 1 matrix. Thus, in a completely explicit notation, the 
above matrix is 16 x 16. 

A few remarks are in order. It is immediately clear that the matrices St are all 
lower triangular, since insertion can never decrease loop number. (Recall that Si, which 
is represented by a diagonal matrix, is left out of the algebra). Triangularity makes these 
matrices easy to exponentiate, since the series will cut off at finite order. This property will 
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be crucial in section [7| Triangular matrices are non- invert ible, which makes sense — typically 
each application of St creates many different graphs. There is no unique inverse under St 
of a generic graph (but see section l6?T|l . Finally, by a quick glance in the appendix we see 
that the matrices are very sparse, which is useful to know for computer storage; the size of 
the matrix is the number of relevant diagrams squared, which quickly becomes prohibitive 
if sparsity is not exploited. 

It will be useful in the following to consider exponentiating matrices, for example the 
matrix representing s(^) in (fT2|) : 

exp(5(^))X = {l + s{yi) + lsi^)smX 

12 2 2 

since the matrix (fT^ vanishes when cubed. On the other hand, 

exp{s{SL))X = . 
Thus, the exponential of the sum of all grafting operators acts as 

where the dots denote higher-order terms. Acting with the inverse is now a simple matter: 

e-S.^*X = X-^X+^X + ^X + ^S + ... (13) 

With the three-loop matrices given in the appendix, one easily performs exponentiation to 
three-loop level as well (see section IH)). The alternating signs in equation (fT^ is already now 
suggestive for the reader who is familiar with Hopf algebra. The antipode of the Hopf algebra 
has a similar alternating sign — each shrinking of a subgraph comes with a sign change. In 
section [7| we display the relation to the antipode, and thus the counterterms of the theory. 
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6 Other Operations 



6.1 Shrinking 



When acting with the transpose of the matrix representation of a grafting operator, we 
shrink graphs. Here is an example: 



/0 1/3 0\/0\ 

2/3 2/3 2/3 0(0 

1 





VOO 0/Vo/ 



/ \ 

2/3 




V / 



= 2/3^ 



This operation does not have a counterpart in the star notation; we will not explore it 
further here. 



6.2 Gluing and External Leg Corrections 

By insertion we only create 1-particle- irreducible (IPI) diagrams. We can also define the 
operation of gluing, that is, joining any two external legs of two graphs with a new propagator. 
Under this operation, either the insertion graph becomes an external leg correction to the 
object graph (if the insertion is has 2 external legs) or we add more external legs (if the 
insertion has n external legs with n > 2). In the latter case, we obtain n-point functions 
with n > 4, which are not superficially divergent in 0^ theory, so we do not consider them. 
As for the external leg correction, it does not have to be considered here either, since we 
only need amputated graphs for the ^-matrix. 

Thus, the gluing operation will not be of much direct interest to us, but it does have 
one property we wish to emphasize. Define ti o to be gluing of ti onto t2, e.g. ^ oiL is an 
external leg correction to ^. The operator s satisfies the Leibniz rule on these diagrams: 

S(tl)(t2 O h) = (S(ti)t2) O t3 + ^2 O isih)ts) . (14) 

It is easy to check that this relation holds in examples in 0^ theory. 

The structure of (jH)) also seems to indicate a possible interpretation as a coderivation. 
We will not explore this direction in this work, since as we have seen, the gluing operation 
is not needed for computing S'-matrix elements in our model field theory. 
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7 Renormalization 



Consider a matrix Mi depending on a parameter e, with elements that diverge when e — > 0. 
Writing Mi(l/e, e), we mean that both positive and negative powers of e occur in the matrix 
elements. If only e occurs in the argument and not 1/e, there are no negative powers. We 
can decompose the divergent matrix Mi into a finite part M2 and a divergent part M3 by 
the general procedure known as Birkhoff decomposition: 



which is uniquely fixed by a boundary condition M2(0) = 0. The matrix M3 has a limit 
as e — s> 0. Birkhoff decomposition was applied to renormalization by Connes and Kreimer 
|2Uj . We take a somewhat different approach; those authors treat the three matrices in the 
decomposition (fT3j) as the same type of objects, whereas we will use the Lie algebra to reduce 
Ml and M3 to vectors, but keep M2 as a matrix. The main point is that we choose M2 in 
the group corresponding to the Lie algebra of matrices defined in section 15.21 That is, the 
divergent matrix M2 is a matrix C of the form 



with numerical coefficients C{t) representing the overall divergence of the graphs t. In 
dimensional regularization, e represents t = A — D where D is the (complex) spacetime 
dimension. The boundary condition M2(0) = is realized in the MS scheme, but other 
schemes can be accommodated by adapting other boundary conditions. The coefficient C{t) 
will depend on e so it should properly be denoted Ce(t), however, we will suppress this 
dependence. The coefficient C{t) is a polynomial in 1/e with no constant term (again, this 
is the boundary condition for ()15|l ). but can depend on external momenta. To calculate the 
complete set of counterterms, it suffices to know these overall- divergence coefficients. 

Let us describe the renormalization procedure in this framework, assuming we know 
the matrix C, i.e. both the overall- divergence values C{t) and the combinatorial matrices St- 
Denote the bare value of a graph t by B{t) and the renormalized graph by A{t) (in the graph- 
by-graph method JE])- Vectors containing these values, indexed by graphs, are denoted A 
and -B, respectively. We Birkhoff-decompose the vector of bare values as in equation (fT3j) : 



Mi(l/e,e) =M2(l/e)M3(e) , 



(15) 
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To find renormalized values from bare values, we liave to invert tlie matrix C, wliicli is a 
trivial matter when it is expressed as an exponential: 



A, 




-Y^C{t)s{t)]By,,, 



(16) 



This is our main statement. As we already pointed out, the sign in the exponential reproduces 
the sign of Zimmermann's forest formula [2], since every s{t) factor in the expansion of the 
exponential corresponds to an insertion, and a sign. 

This is most easily understood in an example. Let us calculate the renormalized 
4-point function up to two loops: 



exp -5^C(t).(t) 





1 














o\ 




(BiX)\ 






1 
















Bm 




-C(X) + 2C(X)C7(X) 




1 













Bm 




-C{^) + C{)^)C{'$} 


-cm 





1 










Bm 




-c(:g[) + 


-cm 








1 







Bm 


v 




-CiSL) 











1/ 




\ Bm J 



where we, as mentioned earlier, it is convenient to use a symmetric renormalization point so 
that C(^) = C()o() = C(!s!). In fact, these three graphs always appear in the same place 
in the expansion, since = s()o() = s(!q!). We should also emphasize that juxtaposition 
of C{t) cannot be interpreted as multiplication. Instead, looking at the matrix C~^, a term 
such as C()o()C(^) corresponds to (l/2)(()o()^) where the angle brackets ( ) pick out the 
pole term in minimal subtraction (in general, the 1/2 is 1/iV!). 

We have, for the nontrivial elements. 



(Am] 






-C(X)5(X) H 




Am 




-{cm - 


-2cmc{yi))B{x) 


-2cmB{-^)+Bm 


Am 




-{cm 


-Ciyx)cm)B{X) 


-cmB{yi)+Bm 


Am 




-{c{^) 


-ci^)cm)B{x) 


-cmB{yi)+Bm) 


\ Am I 




\ -{c{^) 


-cm)c{sL))B{x) 


-C{SL)Bm + B{^) J 



To be specific, consider the third row: 

^(^) = - c{y^)cm)B{x) - c{yi)Bm + • 



(17) 



18 



The combination of graphs is clearly correct; the second counterterm (third term) cancels the 
potentially nonlocal divergence in and the first counterterm (first and second terms) takes 
care of the overall divergence after the nonlocal one is canceled. There is an even number 
of terms, as expected from the underlying Hopf algebra structure. The matrix C "knew" 
the connection between a graph and its divergent subgraphs, since C was constructed from 
the grafting operators St- We can check the cancellation of nonlocal divergences in ()17|) . 
Evaluating now in (p^ theory, with /i, /2, gi, g2 containing no negative powers of e and at 
most polynomial in p: 



2 2 1 

— - -lnp^ + -/i) - (^- - Inp"" + -gi 



2 2, . 2 \ /2 2, . 1 , , 

---In/ + -(^, + 2(72 + ^--ln/ + -/i + /2 

r(/2-2<72). 



2(47r)2" 

Here we included the symmetry factors (see Appendix^ of 1/2. We see that the cancellation 
of (l/e)lnp2 is assured. The second check is that the second order vertex counterterm 
(C(^)-C(W)C(X)) (X (2^?i - /2)/e is local. 

The total contribution to the 4-point function up to two- loop order is just the sum of 
all compatible (4-external-leg) elements of A, as in 

A = A{X) + A(X) + A{X) + A{^) + A{^) (+crossing) 

The total combination of graphs can be read off from the A vector above. 



8 Three-loop Example 

All the diagrams at 3-loop order in (f)^ theory have been computed, and clever techniques 
have been invented to perform the computations. Some techniques are reviewed in ^Hj. A 
useful reference for tables of divergent parts up to three loops is ^Sj. Therefore, we need 
only worry about the combinatorics, but for completeness we try to convey an idea which 
integral computation techniques can be used to compute the relevant diagrams. 

In one-loop computations, one usually employs the method of Feynman parameters. 
In fact, using Feynman parameters without additional techniques is, by some estimates, 
only economical for graphs with two lines. "Additional techniques" can include performing 
some series expansions before computing the integral (the Gegenbauer polynomials have 
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been shown to be well suited for this), and imaginative use of integration by parts. In this 
context, integration by parts is applied to simplify integrals using the vanishing of integrals 
over total derivatives: 

for any integrand I{p) depending on the loop momentum p and any other loop or external 
momenta. Applying integration by parts, the massive ladder diagram ^ at zero momentum 
transfer can be decomposed as jjl] 



^ (47r)4 V47rM2 J J 



1 

+ - 

e 



(27r)^ (F + m2 

(47r)2 \4nM^ ) 



where F is a simpler genuine three-loop integral, i.e. it cannot be expressed in terms of 
integrals of lower loop order, Wq is a two-loop integral at zero momentum transfer) and 
5*3 is a one- loop integral at zero momentum transfer). The result (given in ^5]) is 



4 / 2 \ -3e/2 

m I m \ 



(An) 



47 

3 ' 



+- 



TT 



--4a -47 + 37" + 
3 o 



+ (finite) 



where 7 is Euler's constant, a is a numerical constant (a = 1.17...) coming from integration 
over Feynman parameters, and M is the renormalization scale. 

Now let us use the matrix Lie algebra to compute the renormalized graph (again work- 
ing in the graph- by-graph context). From the previous section, we know that counterterms 
are generated by the exponential of a sum over grafting operators. 

Using 



A = exp 



(c^(X)s(X) + c{x)s{X) + c^(^)s(^) + ■■■] 



B 



we find the renormalized graphs. In particular, the relevant row of the vector A is 

^(M) = (-^C'(X)' + ^C(X)C(^) + ^C(X)C(X)-iC()§())i?(X) + 
iC^(X)'5(X) - |C^(X)5(S) - \c{:$)B{X) + bC^) . 
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This is to be compared to the known result for the renormahzed graph: 

-()o()X-(>o()^+()o()2)o( (19) 

where ( ) denotes the renormalization map (for example, in minimal subtraction we simply 
drop the finite part). Rewriting the known expression (|TU|) in language more similar to the 
above, we have 

Ai^) = i?()§() + (-C()§() + C(X)C(^) + C(X)C(X)-C(X)^C(W))i?(X) 
-C()o()E(^) - C{)oi)B{-^) + C()o()25()o() . 

Expression ()19p is reproduced by that derived by the Lie algebra, up to the differing normal- 
ization (see section E}. To summarize, we have used the inverse exponential of Lie algebra 
elements to generate counterterms, just as the antipode would have been used in a Hopf 
algebra. 

The Ihara bracket relation siti,t2] = [■^ti, s^j] is still rather trivial in this example, as it 
turns out, because only s(^) appears more than once. In the example of we have s(^), 
s(X) and sC^) appearing, so there is potential for use of the relation here. However, the 
only nontrivial commutator, equation (@)), yields graphs that do not appear as subgraphs of 

i s(^) + ^ si^) - 2 s(^) = [s(X), sin)] . 

We have thus seen that it is not until the four-loop level that nontrivial application of the 
Lie algebra relation S[tj^^t2] = [^h, Stj] appear. As mentioned in the introduction, new results 
are needed at five loops and higher. 



9 Renormalization Group Flows and Nonrenormaliz- 
able Theories 

In Wilson's approach to renormalization, short-distance (ultraviolet) fluctuations are inte- 
grated out of the functional integral to appear only as modifications of parameters of the 
Lagrangian for the remaining long-distance (infrared) degrees of freedom i37j . The introduc- 
tion of a mass scale M by renormalization parameterizes different renormalization schemes, 
and the change of M induces a flow of parameters in the Lagrangian, such that certain 
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so-called irrelevant operators die away. In 0^ theory in four dimensions, 0^ is an example of 
one such operator; the relative size of this operator to other terms in the Lagrangian at a 
momentum scale p would be {p/ K)^"^ = (p/A)^ as p — > 0, where A is a cutoff. 

Now, the grafting operator St may be thought of as a "magnifying glass" ; by inserting 
graphs, we resolve details of graphs that were not visible at larger scales. In this specific 
sense, St induces scale change. In particular, in this paper, we only considered Sj for graphs t 
with 2 or 4 external legs. We do not have to consider st for graphs with 6 external legs in (f)'^ 
theory in four dimensions, precisely because 0® is an irrelevant operator in this theory. This 
shows how the matrix Lie algebra would appear in a nonrenormalizable theory; everything 
is the same as in the 0^ example, except there is an infinite number of insertion matrices. 
While this situation is certainly more cumbersome, it is not necessarily fatal. 

According to Connes and Kreimer (20], renormalization group flow boils down to the 
calculation of a matrix (3 which in our notation becomes 

t 

and is independent of e. Here I3{t) is the numerical contribution to the beta function due to 
a certain diagram t (again, in the graph-by-graph method). This matrix (3 generalizes the 
standard beta function of the renormalization group; there is now combinatorial information 
attached to each of the contributions I3{t) to the beta function. 

See also the next section for some comments on the range of validity of Wilson's point 
of view. 

10 Conclusion and Outlook 

In this work, we displayed a matrix Lie algebra of operators acting on Feynman graphs, 
that by exponentiation yields group elements generating counterterms for diagrams with 
sub divergences. We defined a grafting operator st that inserts one Feynman graph t into 
another. The matrix representations of these operators satisfy a certain rule, similar to 
an Ihara bracket, that gives relations between them: s\t^^t2\ = ["^tuS^j]. In this way, not all 
matrices have to be separately computed, with potentially substantial savings in computation 
time. We displayed the relation to the star product of Kreimer, which is defined similarly 
to (but not exactly the same way as) the grafting operator. A simple computation verifies 
that the right-symmetric (pre-Lie) nonassociative algebra of the star product is equivalent 
to the previously mentioned Ihara-bracket rule for our matrix representations. 
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We also gave a number of examples, mostly rather simple ones, and checked in a three- 
loop example that the correct sequence of counterterms is provided by the exponential of 
Lie algebra elements. (The general proof that this is always the case will be given elsewhere 
inj). Just as with Hopf algebra, the Lie algebra rules are trivial at one-loop, almost trivial at 
two-loop, and just beginning to become interesting at three-loop order. At four loops, there 
is plenty of interesting structure to check, but it will require a fair amount of computation; 
this is an obvious direction in which future work should go. 

An equally obvious, but less direct, application is to noncommutative geometry. Using 
the Connes-Kreimer homomorphism between the Hopf algebra of renormalization and the 
Hopf algebra of coordinates on the group of formal diffeomorphisms pUl , and the Milnor- 
Moore theorem, the Lie algebra in this paper has a corresponding Lie algebra in the context 
of noncommutative geometry. We hope that the matrix Lie algebra may eventually shed 
some light on certain aspects of noncommutative geometry, in this rather indirect but novel 
way. 

In a less obvious direction, it is suggestive to note that the grafting operator St is an 
element of a Lie algebra, it satisfies a Leibniz rule (the gluing operator, section IHS]), and its 
bracket (jH)) looks like that of a vector field X on a manifold. Loosely, the grafting operator is 
a "scale vector" , that takes us into ever increasing magnification of graphs by including more 
and more subgraphs. If a Lie derivative Cgt along this "vector field" could meaningfully be 
defined, that would open the following interesting speculative possibility. 

It was proposed in 0, based on earlier work [7^ that volume forms for functional 
integration be characterized the following way: find a vector field X on the space, define 
its divergence (this may be a nontrivial task in an infinite-dimensional space) and let the 
volume form uj be defined by 

Cx = (divX) uj . 

This definition reproduces familiar volume forms in some simple finite-dimensional examples 
(Riemannian manifolds, symplectic manifolds), and gives some hope for generalization to 
infinite-dimensional spaces. Perhaps if St is a vector field, it could be used in the role of X, 
in some extended sense, to characterize u. In effect, this would be a perturbative definition, 
since including St up to a certain loop order will define different u at different orders in the 
loop expansion. 

An obvious problem with this idea is that St acts on the space of graphs, not the space 
of fields. This means that even if a volume form u could be defined, it would be a function 
on the space of graphs. In principle, it may be possible to exploit some analogy to the space 
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of paths in quantum-mechanical functional integrals. This direction would be interesting to 
pursue in future work. 

As a final note, there has been some dispute to what extent Wilson's picture is gen- 
erally valid; there are some (as yet speculative) examples where infrared and ultraviolet 
divergences are connected jBHI- Some of these examples are in connection with noncommu- 
tative geometry. In view of the relation between the algebraic approach to renormalization 
and noncommutative geometry mentioned above, one could hope that these algebraic devel- 
opments may eventually shed some light on the ultraviolet-infrared connection. 
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A Symmetry factors 



Since conventions for how graphs carry symmetry factors vary, let us be precise. 

The bare graph in our two-loop example contains in its bare value the nonlocal diver- 
gence (l/e)lnp^, including the 1/2 due to the symmetry factor 2. The counterterms must 
subtract the (1/e) Inp"^, nothing more, nothing less. 

How the notation carries symmetry factors, however, varies in the literature. The 
term ()o( )^ taken literally looks like it would contain two symmetry factors of 2, i.e. the 
forefactor would be 1/2 x 1/2 = 1/4. By itself, this is of course wrong — it only subtracts 
half the nonlocal divergence. Two possible resolutions are i) to say that the graph with a 
cross (as in JE]) that represents ()o()^ has only symmetry factor 2, not 4 (which is not quite 
consistent notation on the non-graph level without a factor of 2, if ( t ) is to really mean 
the pole part of t), or ii) to posit, as we do here, that insertion can also create a term with 
C(5s^), that combines with C()o() when we take the counterterms to be the same, to give a 
factor of 2. 

It is a useful check to consider the s-channel counterterm as well. In contrast to the 
previous example, the graph with an s-channel insertion (the first L = 2, E = 4 graph in 
our table 1) has a forefactor of 1/4. Each counterterm also has 1/4, but now there are 
two of them, one for shrinking each of the two disjoint sub divergences. Thus they do not 
cancel directly, but factor as usual. In the Hopf-algebra coproduct, factorization produces 
the factor of 2. The antipode then yields a factor of 2 similar to the above. With these 
conventions, finding a factor of 2 both in the s and t/u cases, one would expect the s/t/u 
graphs to have the same normalization, after summing. 

B Insertion Tables 

In this appendix, we provide tables of the action of the grafting operator St- We have 
defined St such that on a bare diagram denoted by 1, we have Stl = t. The matrix si acts 
as Sit = N{t)t where N{t) is the number of insertion points of t. It is also true that the 
coefficients of each resulting graph in an insertion into a graph t should add up to the number 
of compatible insertion points of the graph t. 

One-loop: 
1 + = 1 
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S{IL)- 

Two-loop: 
1 + 1 = 2 

s{SL)SL 

Three-loop: 
1 + 2 = 3 

s(iL)_8_ 






- + crossing 



_8_ 



2^ 

1 2 

2 2 2 

3X+3S+3K 



2^ 

2 \^ X 2 2 v->v 2 

-$+-$+-5:+3g+3g 
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2 + 1 = 3 



s(_a_)iL = X 

s{^)^ = 2^ 
s{X)SL = ^±+'^^ 



C Grafting Matrices 



In this appendix, we list matrix representations of the grafting operators s{t). Since insertion 
into an n-point function can only create an n-point function, we consider submatrices of 
each separately, and denote them by subscripts. For example, if we call s(X) = A and 

— B, then the complete matrix representing s(^) is the direct sum of A and B: 




Lines are drawn through the matrices to delineate loop order L. 
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One-loop insertions: 
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When it comes to three-loop matrices, it is typographically easier to give the trans- 



pose: 



X 



/O 00000000000000000000 o\ 
0000000000000000000000 



2 



10 2 



0000000000000 
OOiiOfOiOfOIOOOOOOOOOO 
00ii0020i0202000000000 

\ooooooooooooo2ioo2||ooy 



X 
X 



where all elements below those displayed are zero (it is, of course, a square matrix), and 
we have not displayed the zero-one-two-loop submatrices already given. That is, the given 
matrix is the transpose of the 22 x 6 submatrix B in 



A 





B 






where A is the 6x6 matrix given for s(^) earlier. The order is as in tabled] summarized 
here: 
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We also give two two-loop insertions: 



X 



/O 00000000000000000000 o\ 
0000§§§000000000000000 
0000000000000000000000 
0000000000000000000000 
0000000000000000000000 

\ooooooooooooooooooooooy 



X 
X 



X 



/O 00000000000000000000 o\ 
OOOiOOOiiiiiiOOOOOOOOO 

3 6 6 3 3 3 3 

0000000000000000000000 
0000000000000000000000 
0000000000000000000000 
\0000000000000000000000/ 



X 
X 



We note that similarly to the previous case, there is also a top left 6x6 matrix with entries 
in the first column only. The rest of the matrices are now trivial to extract from the insertion 
tables, so we shall not repeat them here. 



References 

[1] D. Becirevic, Ph. Boucaud, J. P. Leroy, J. Micheli, O. Pene, J. Rodn'guez-Quintero, 
C. Roiesnel, Asymptotic scaling of the gluon propagator on the lattice, Phys.Rev. D61 
(2000) 114508; hep-ph/99 10204| 

[2] N.N. Bogoliubov, O.S. Parasiuk, Uher die Multiplikation der Kausalfunktionen in der 
Quantentheorie der Felder, Acta Math. 97 (1957) 227; K. Hepp, Proof of the Bogoliubov- 
Parasiuk Theorem of Renormalization, Commun. Math. Phys. 2 (1966) 301; W. Zim- 
mermann, in S. Deser et. al. eds. Lectures on Elementary Particles and Quantum Field 
Theory, Vol I (MIT Press, Cambridge, 1970). 

[3] D.J. Broadhurst, D. Kreimer, Renormalization automated by Hopf algebra, J. Symb. 
Comput. 27 (1999) 581; |hep-th/98 10087. 



29 



D.J. Broadhurst, D. Kreimer, Combinatoric explosion of renormalization tamed by 
Hopf algebra: 30-loop Pade-Borel resummation, Phys. Lett. B475 (2000) 63-70; 
fhep-th/QQ 12093, 

D.J. Broadhurst, D. Kreimer, Exact solutions of Dyson- Schwinger equations for iterated 
one-loop integrals and propagator- coupling duality, Nucl. Phys. B600 (2001) 403-422. 

P. Cartier, to be pubhshed. 

P. Cartier, C. DeWitt-Morette, A new perspective on functional integration, J. Math. 
Phys 36 (1995) 2237-2312. P. Cartier, C. DeWitt-Morette, Functional integration, J. 
Math. Phys 41 (2000) 4154-4187. 

P. Cartier, M. Berg, C. DeWitt-Morette, and A. Wurm, Characterizing Volume Forms, 
to appear in Hagen Kleinert's festschrift (World Scientific, 2001); |math-ph/0012009, 

A.H. Chamseddine, A. Connes, The Spectral Action Principle, Commun. Math. Phys. 
186 (1997) 731-750. 

F. Chapoton, Un theoreme de Cartier-Milnor-Moore-Quillen pour les bigebres dendri- 
formes et les algebres braces; math.QA/0005253,. 

K. Chetyrkin, M. Misiak and M. Miinz, Beta functions and anomalous dimensions up 
to three loops, NucL Phys. B518, 473 (1998); |h^ph/971 1266 



K.G. Chetyrkin, T. Seidensticker, Two Loop QCD Vertices and Three Loop MOM j3 
functions, Phys.Lett. B495 (2000) 74-80; pIep^/0008094 

K.G. Chetyrkin, F.V. Tkachov, Integration by parts: the algorithm to calculate beta 
functions in 4 loops. Nucl. Phys. B192 (1981) 159-204; K.G. Chetyrkin, A.H. Hoang, 
J.H. Kuhn, M. Steinhauser and T. Teubner, QCD corrections to the e^e~ cross section 
and the Z boson decay rate: concepts and results, Phys. Rep. 277 (1996), 189-281. 



C. Chryssomalakos, H. Quevedo, M. Rosenbaum, J.D. Vergara, |hep-th/0105259 



J.-M. Chung and B.K. Chung, Calculation of a Class of Three-Loop Vacuum Diagrams 
with Two Different Mass Values, Phys. Rev. D59 (1999) 105014; Tiep-ph/9805432^ 

[16] J. C. Collins, Renormalization (Cambridge, 1985). 



30 



[17] A. Connes, A Short Survey of Noncommutative Geometry, J. Math. Phys. (2000) 
3832-3866; hep-th/0003006, 

[18] A. Connes, M.R. Douglas, A. Schwarz, Noncommutative Geometry and Matrix Theory: 
Compactification on Ton, J. High Energy Phys. 02 (1998) 003; |hep-th/9711162[ 

[19] A. Connes, D. Kreimer, Renormalization in quantum field theory and the Riemann- 
Hilbert problem I: the Hopf algebra structure of graphs and the main theorem, Commun. 
Math. Phys. 210 (2000) 249-273; ,hep-th/9912092, 

[20] A. Connes, D. Kreimer, Renormalization in quantum field theory and the Riemann- 
Hilbert problem II: the f3-function, diffeomorphisms and the renormalization group, 
Commun. Math. Phys. 216 (2001) 215-241; |hep-th/0003188| 



[21] A. Connes, D. Kreimer, Lessons from Quantum Field Theory — Hopf Algebras and 
Spacetzme Geometries, Lett. Math. Phys. 48 (1999) 85-96; hep-th/9904044j 

[22] A. Connes, H. Moscovici, Differentiable cyclic cohomology and Hopf algebraic structures 
in transverse geometry; math.DG/0102167| 

[23] J.M. Gracia-Bondi'a, J.C. Varilly, H. Figueroa, Elements of Noncommutative Geometry 
(Birkhauser, Boston, 2001). 

[24] T. Hahn, Generating Feynman Diagrams and Amplitudes with FeynArts 3, 
|hep-ph/0012260| The FeynArts web site is www . f eynarts . de 



[25] S. Heinemeyer, Two-loop Calculations in the MSSM with FeynArts, hep-ph/0102318,. 
[26] C. Kassel, Quantum Groups (Springer- Verlag, New York, 1995). 

[27] M. Kontsevich, Operads and Motives in Deformation Quantization, Lett. Math. Phys. 48 
(1999) 35-72; math.QA/9904055l 

[28] D. Kreimer, On the Hopf algebra structure of perturbative quantum field theories. Adv. 
Theor. Math. Phys. 2 (1998) 303-334; q-alg/9707029j 

[29] D. Kreimer, Shuffling Quantum Field Theory, Lett. Math. Phys. 51 (2000) 179-191; 
|hep-th/99 12290] 

[30] D. Kreimer, Combinatorics of (perturbative) quantum field theory, to appear in Phys. 



Rep.; hep-th/00 10059 



31 



[31] D. Kreimer, Knots and Feynman Diagrams, (Cambridge, 2000). 



CP. Martm, J.M. Gracia-Bondia and J.C Varilly, The Standard Model as a noncom- 
mutative geometry: the low energy regime, Phys. Rep. 294 (1998) 363-406. 

[33] S. Minwalla, M. Van Raamsdonk, N. Seiberg, N oncommutative Perturhative Dynamics, 
J. High Energy Phys. 02 (2000) 020; hep-th/9912072; M. Van Raamsdonk, N. Seiberg, 
Comments on N oncommutative Perturhative Dynamics, J. High Energy Phys. 03 (2000) 
035; hep-th/0002186; L. Griguolo, M. Pietroni, Wilsonian Renormalization Group and 
the Non- Commutative IR/UV Connection, ,hep-th/0104217|[ 

[34] T. Ohl, Drawing Feynman Diagrams with LaTeX and Metafont, Comput. Phys. Com- 
mun. 90 (1995) 340; |hep-p h/9505351| 

[35] G. Racinet, Series generatrices non-commutatives de polyzetas et associateurs de Drin- 
feld, Ph.D thesis (Universite de Amiens, France, 2000). 

[36] N. Seiberg, E. Witten, String Theory and N oncommutative Geometry, J. High Energy 
Phys. 9909 (1999) 032; |hep-th /9908142| 

[37] K.G. Wilson, J. Kogut, The Renormalization Group and the e Expansion, Phys. Rep 
12 (1974) 75-200. 



32 



